Stable Skyrmions in SU{2) Gauged Bose-Einstein Condensates 
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We demonstrate that the three-dimensional Skyrmion, which has remained elusive so far, spon- 
taneously appears as the ground state of SU (2) symmetric Bose-Einstein condensates coupled with 
a non-Abelian gauge field. The gauge field is a three-dimensional analogue of the Rashba spin-orbit 
coupling. Upon squashing the SO(3) symmetric gauge field to one- or two-dimensional shapes, we 
find that the ground state continuously undergoes a change from a three-dimensional to a one- or 
two-dimensional Skyrmion, which is identified by estimating winding numbers and helicity. All of 
the emerged Skyrmions are physically understandable with the concept of the helical modulation 
in a unified way. These topological objects might potentially be realizable in two-component BECs 
experimentally. 
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Introduction. — The n-dimcnsional Skyrmions (n<3), 
classified by the n-th homotopy group (ir n (S n ) =Z), have 
attracted much attention in various research fields, rang- 
ing from high-energy to condensed-matter physics 13|] - 
It has been demonstrated that the two-dimensional (2D) 
Skyrmion spontaneously appears as the ground state in 
helical magnets MnSi and Fei-zCo^Si [4|-|6|], a quan- 
tum Hall state 0, and continuous vortices analogous to 
Skyrmions in 3 He- A Furthermore, it has recently 

been created by using the phase-imprinting technique in 
gaseous Bose-Einstein condensates (BECs) [l(|. So far 
those are all 2D objects. 

The three-dimensional (3D) Skyrmion is a particlelike 
soliton in classical field theory, which was hypothetically 
introduced by Skyrme to describe baryons in a meson 
field. Although this hypothesis has attracted a lot of at- 
tentions for decades, the evidence has yet to be clarified. 
One difficulty of the proposal is due to the instability of 
Skyrmions: it is known that in the nonlinear cr-model, 
the gradient energy makes the 3D Skyrmion unstable to- 
ward shrinkage, in spite of the topological stability Q 
because energetics differs from topology. To prevent it 
from shrinkage, Skyrme 0] added by hand the quartic 
differential term, which is the so-called skyrme term. In- 
stead of adding an ad hoc Skyrme term, it has been clar- 
ified that a non-Abelian gauge field in the form of the 
t'Hooft-Polyakov monopole, which yields the scaling law 
the same as the skyrme term, facilitates the stability of 



the 3D Skyrmion It is, however, a nonrealistic and 



purely theoretical proposal, because such a non-Abelian 
gauge field does not exist in the meson field theory 

On the other hand, the stability of the 3D Skyrmions in 
multicomponent BECs has been investigated since their 
proposals II, 13l. Although the schemes to create and 
stabilize them |l44ig ] have theoretically been proposed, 
the 3D Skyrmion is still elusive both experimentally and 
theoretically. One reason is that the Skyrmions in previ- 
ous works are merely metastable solutions of the energy 



functional 14- 19j|. 

The aim in this Letter is to clarify that a 3D Skyrmion 
spontaneously emerges as the "ground state" of BECs, 
coupled with a realistic non-Abelian gauge field, with- 
out the help of the Skyrme term. This is the first pro- 
posal of the stable Skyrmions with the 3D analogue of 
Rashba SOCs. First, we demonstrate that the stability 
of the 3D Skyrmion is physically understandable with 
the concept of the helical modulation of the order pa- 
rameter (OP) [13, [2l|| . We show the phase diagram and 
the stable Skyrmion textures by numerically solving the 
full Gross-Pitaevskii (GP) equation. Here, we mainly 
focus on a BEC with an SU(2) symmetric interaction 
to capture the essential physics in the presence of the 
non-Abeli an gaug e field, as widely studied in the ear- 
lier works [l44l9| . The Hamiltonian is analogous to the 
Higgs sector of the Weinberg-Salam model of electroweak 
interactions (22|. However, we also consider the stabil- 
ity of the 3D Skyrmion against the interaction without 
SU(2) symmetry. 

Recently, there was a major breakthrough that en- 
abled one to artificially imprint a gauge field in ultracold 
atoms [23l - [2"5l | . This is intriguing in the sense of acces- 
sibility to new topological phases and the appearance of 
spatially modulated ground states due to non-Abelian 
gauge fields j2^|. The technique is based on the Raman 
coupling between hyperfinc states of atoms which recon- 
structs the internal degrees of freedom to be degenerate 
pseudospin states. The Hamiltonian in the pseudospin 
representation has a fictitious Abelian or non-Abelian 
gauge field due to the adiabatic motion of its degener- 
ate pseudospin states. 

The schemes to generate various types of gauge fields 
have been proposed theoretically. In two-component 
BECs, there exist methods to generate 2D and 3D ana- 
logues to the Rashba spin-orbit coupling [29l - [3~D ] and a 
monopole field with a Dirac string [32j. Recently, the 
one-dimensional (ID) Rashba+Dresselhause type gauge 
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field was experimentally realized by the NIST group [25| . 
Fermi gases coupled with a non-Abelian gauge field were 
also investigated [2(| [2?]]. In this Letter, we clarify how 
such non-Abelian gauge fields stabilize the 3D Skyrmion. 

Phase diagram. — We start with the GP energy func- 
tional for two-component bosons with the atom mass m 
and pseudospins f, 4- as 



E = H 
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where the particle density is n(r) = and the z- 

componcnt of the local spin is S z (r) = ■^ 1 i'*^(cr z ) f j jl , 1 if 1/ . 
Throughout this Letter, the repeated Greek indices im- 
ply the sum on the spin /x, v, r/ =t, I and we use the 
unit of h = m = u) = 1, where oj is the trap frequency. 
The Co + ci is the interaction strength between bosons 
in the same pseudospin component and the Ci indicates 
the interaction between the intercomponents. The OP 
<J M for the condensate obeys the normalization condi- 
tion, J 'I'* iff^cPr = 1. The single particle Hamiltonian 
Hq in Eq. (JTJ is given by 

Ho = J d 3 r [D^Hfuirtf ■ [A*,*„(r)] . (2) 

Here, we set the covariant derivative -D^y = 
(— iVcro + -A) with the 2x2 unit matrix oq and a non- 
Abelian gauge field A. In most of this Letter, we set 
ci = in Eq. (jT]) . which ensures the SU{2) symmetry of 
the interaction. However, we will mention later that the 
finite ci does not make the Skyrmion unstable. 

First of all, we summarize in Fig. QJa) the schematic 
phase diagram of BECs coupled with the 3D non-Abelian 
gauge field, 



A = kj_ (a x x + a y y) + K z a z z, 



(3) 



where <jj denotes the j-th Pauli matrix. This is 3D ana- 
logue of a Rashba or Drcssclhause type spin-orbit cou- 
pling which is known in the condensed matter context. 
In fact, the scheme to generate this type of synthetic 



gauge field with k z = k± [3l| and k z = [29|, |30( was 
proposed theoretically, and the case of k± = was real- 
ized in the experiment [25| as mentioned. In Fig. QJa) , 
we find out the stable region of the 3D Skyrmion tex- 
ture, while the other regions are occupied by the ID or 
2D Skyrmion. It turns out that these textures smoothly 
change to each other. It is demonstrated in the rest of 
this Letter that all types of Skyrmions, whose textures 
are displayed in Fig. [2] are understandable with the he- 
lical modulation of the OP which can be parameterized 
on the three-dimensional surface S 3 in the OP space. 

Helical spin modulation and Skyrmions. — The concept 
of the helical modulation of the OP 21( provides a good 
starting point to the understanding of the stability and 
smooth transition of Skyrmions under the non-Abelian 
gauge field described in Eq. ©. 
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FIG. 1: (Color online) (a) Schematic phase diagram of tex- 
tures and the symmetries of Hq. At n z — n±, the Hamiltonian 
has an SO^n+s symmetry, (b) Schematic picture of the heli- 
cal spin modulation along h. The thick green (gray) arrow de- 
picts the local spin S\ 0) = ±{¥°\r )}^a ] ) tl M°' ) (r )/n(r ). 



The OP of a two-component BEC is always parameter- 
ized, the U(l) phase <I> = $(r), and the angles s = (/> s (r) 
and 8 s = 6 s (r), as 



* t (r) 



-i<Tz B /2 -i<J v 6 B /2 



(4) 



Here, as shown in Fig. (Ub), <j> s (r) and 8 s (r) de- 
note the direction of the local spin, Sj(r) = 
i** {r)(aj)^ v {r)/n(r). Since the OP manifold in 
Eq. ((3]) is mapped onto a three-dimensional surface 
SU(2) ~ S 3 , two-component BECs have a topological 
object classified by the homotopy group 7r 3 (S' 3 )=Z. 

Turning now to a non-Abelian gauged BEC, we demon- 
strate that the helical modulation of the OP due to the 
field A makes the Skyrmion stable even without the 
Skyrme term. It is crucial to observe that the helical 
modulation with an arbitrary modulation vector h can be 
written down within the SU (2) symmetric OP in Eq. (|4]) 
as 



^(r,h)=U^(h,2h-r)^\r ), 



(5) 



where the 2x2 matrix U{n, (p)=exp[i((p/2)(n-(r)] £SU(2) 
denotes the SU(2) rotation of an arbitrary OP 
around n by the angle tp. Figure [ljb) shows the 
schematic picture of this modulation. The helicity origi- 
nates from the single particle spectrum Eq — Hq of ideal 
Bose gases in the thermodynamic limit Vl/^ (r ) = e lk r ^^ , 
which is given by E a = k 2 + n 2 ± 2^J K 2 z k1 + + k%) 



with 



■2k\ 



It turns out that the OP of the 



ground state is spatially modulated, since Eq has min- 
ima on the finite fc. Then, the helical modulation with h 
in Eq. ([5]) is the superposition of momentum eigenstatcs 
k || ±h. The spatial inversion symmetry of Jio guarantees 
the degeneracy of the helical modulation starting from an 
arbitrary Vl/^ . 

For k z = k± = n/y/3, the single particle Hamiltonian 
Hq in Eq. ([I} is invariant under the simultaneous ro- 
tation of spin and real spaces 5*0(3)1^48. Since Eq has 



3 



minima on surface k = /t/Vo, the helical spin modula- 
tion in Eq. ([5]) is degenerate for any direction of h\\ k. 
Note that the 3D helical modulation ^^(r,h |j r) prop- 
agating with all the directions of h along r fulfills the 
OP manifold S 3 within < h ■ r < tt. Thus, this texture 
v]/3D = j 2h ■ r) is the 3D Skyrmion, which is the 
candidate of the ground state. 

In contrast, for the region k z /kj_ < 1 in Fig. HJa), 
the .S'O^R-fs symmetry in Jio in Eq. (JTJ) is broken into 
the S , 0(2)r4s that denotes the joint rotation of spin and 
real spaces around the £-axis. Eq also has a minimum 
line along k±= k\ + k^ = n±_ . Therefore, it turns out 

that the possible stable texture for k z < k± is the radial 
or the ID helical spin modulation expressed by Eq. ([5]) 
with h || (k x , ky,0), where the former corresponds to a 2D 
Skyrmion [34j . 

For the region k z /k± > 1 in Fig. HJa), the Hamilto- 
nian ([T]) still remains invariant under the discrete sym- 
metry that is the simultaneous 7r-rotation of spin and real 
spaces around the x- or y-axis, where Eg has minima at 
k = ±k z z. Then the most stable modulation vector h 
is confined to the z-axis, implying that the possible sta- 
ble texture is the ID helical modulation along £-axis. 
Since this spin-rotation along the z-axis consists of the 
U(l) degrees of freedom, we can interpret this helical spin 
modulation MS (X ID Skyrmion. 

Stable Skyrmion textures. — In order to quantitatively 
discuss the candidates of the stable spin textures, we 
numerically minimize the full Gross-Pitaevskii energy 
functional ([1]) by using the imaginary time evolution 
scheme [33| with a spatial grid of 121 3 . We have veri- 
fied that the obtained solution is the true ground state 
by starting the calculation from a variety of initial condi- 
tions such as the uniform spin, helical modulation with a 
variety of modulation vectors, and 2D and 3D Skyrmions. 
For numerical calculations, we use Co = 100, where the 
density profile yields a Gaussian-like shape. This pa- 
rameter corresponds to the trap frequency aj/2ir = 100 
Hz and s-wave scattering length of the same component 
ttii = «22 = 5 nm with the 87 Rb mass and number of 
trapped atoms N — 5 x 10 4 . These values define the 
length of the trap unit as y/ H/mui ~ 1 /im. To quan- 
tify the size of the atomic cloud, we use the Thomas- 
Fermi radius Rq = \/2fiQ = 2.99 in the trap unit, where 
/.to = (15co/27r) 2 / 5 /2 is the chemical potential with the 
Thomas-Fermi approximation. 

Figure [2] shows the stable 3D Skyrmion obtained from 
the numerical minimization of Eq. (Q]) in the presence 
of the gauge field with ki = k z and the spherical trap 
potential. Figure. [2ja) shows that the pseudospin texture 
helically modulates along the 3D radial direction. In this 
texture, the singularity of the U(l) phase $ and the angle 
4> s introduced in Eq. (j4j exists on the gray line where 
the local spin points to z (^-singularity) or — z {^^~ 
singularity) direction. As shown in Fig. [^b) and[5Jc), 



(a) ^-ir/2H^n H-jt/2 





FIG. 2: (Color online) The spatial profile of the stable 
3D Skyrmion in the x-y and z-x planes. The arrows and 
their colors in (a) indicate the pseudospin direction and the 
U{\) phase of the OP, respectively. The gray lines in (a) 
imply the singularity of $ and 4> s ■ The color maps of (b) 
and (c) give the amplitudes |$-|-(r)| and |^4.(r)|, respectively. 
The black arrows in (b) and (c) denote the directions of the 
phase winding of each component. These results are ob- 
tained with the parameter (k±Rq, k, z Ro) = (5.07, 5.07), where 
iio = (15co/27r) 1/5 . 

OP component ^ M accumulates the phase 2-7T on the path 
enclosing its singularity. Namely, forms the vortex 
ring and its ring-singularity is fulfilled by Vf^ with the 
phase winding. This can be interpreted as the so-called 
"vorton" structure known in high-energy physics [H l35j . 

The 3D Skyrmion is identified by the winding number 
tt^{S z ) of the map from real space to the OP manifold 
SU(2) ~ S 3 . The winding number defined by Refs. [10- 

[13 

Wm = i / ^reijksmesidiesXdrfsXdk®) (6) 

counts how many times the map warps the OP man- 
ifold. Figure [3][a) shows the winding number W30 in 
the plane of k z /k± and Rok_l, obtained from the nu- 
merical solution of the full GP equation. Here, we esti- 
mate Eq. in the region of r < 1.5Rq. It is seen from 
Fig. [3ja) that Wsd increases with growth of the n±/R 
near the k z /k± = 1 line. Hence, in this region the 3D 
Skyrmion becomes stable. This continuous increase of 
W30 is because of the absence of the boundary condi- 
tion that ^ i_i(r — > 00) is nonzero and uniform, assumed 
for Skyrmions in other contexts. So far as the bound- 



4 




FIG. 3: (Color online) Stereographic plots of W30 (a), 
H20/R' (b), and W\-d/^R' 2 (c) in the plane of k z /k± and 
k±Ro. These denote the winding numbers and helicity which 
quantify the ID, 2D, and 3D Skyrmions. We estimate these 
values in the region of r < 1.5f?o = R' ■ 



ary condition is satisfied, W3D € Z and 3D Skyrmion is 
topologically stable. However, the presence of the vec- 
tor potential A makes ^(r — > 00) nonuniform even if 
we ignore the effect of the trap. In fact, with increasing 
of n z = KA^Rf) 1 , the characteristic length of the helical 
spin modulation becomes smaller and the shell-shaped 
Skyrmion penetrates from outside of the system in addi- 
tion to a unit of the Skyrmion. Therefore, as k increases, 
W3B increases continuously and monotonically as shown 
in Fig. Ufa). 

One- and two-dimensional Skyrmions. — The 2D 
Skyrmion texture is described as VP^ = ^f^(r,h || p), 
where p = x + y. The Skyrmion \E'^ D consists of the 
helical spin modulation along the 2D radial direction, 
which is favored by the gauge field described in Eq. ^ 
as mentioned above. This texture occupies the region 
k 2 /k_l 1 in Fig. QJa). If we ignore the U(l) phase 
the winding number of this 2D Skyrmion can be defined 
as W 2B = ±f d^re l3 sin^c^s)^) Q. 

We should notice that the winding number of vf^ 
is determined by the domain of the rotating angle J7; 
W2D = +1 (— 1) is accumulated within nn < Q, < 
(n+ l/2)7r ((n — 1/2)tt < < mr). This implies that 
with increasing k±, the Skyrmion (W2D=+1) and anti- 
Skyrmion (W2B = —1) penetrates from outside alter- 
nately. Hence, W2D of 2D Skyrmion "J™ oscillates with 
increasing k±. Instead of W20, one can estimate the 
size of the 2D Skyrmion with the helicity introduced as 
H-2u = ^J d 3 r laj sin# s (a^s)(^s)|. 

The ID Skyrmion along the 2-axis is described as 
^jP = *p( r i^ il z )- This appears as the ground state 
in the region k z /k± > 1 of Fig. [TJa), where the modula- 
tion vector h(x±K z z is favored as discussed above. The 



winding number for the ID Skyrmion is introduced as 
Win = 57 J ^z — i£ -^^-^ + c.c. This corresponds to 
the phase accumulation which quantifies the spin current 
of S z along f-axis. 

We plot Win and H 2 d in Figs. [3fb) and (c), where 
they reveal the stable region of the ID and 2D Skyrmions 
in the parameter space spanned by k z and K±. In the 
region k z /k± < 1, the 2D Skyrmion is stabilized where the 
helicity FL 2 v> increases and W3D and Win are suppressed. 
In the other region, k z /k± > 1, the growth of Wm and 
the suppression of W3D and i?2D indicate the appearance 
of the ID Skyrmion. Figure [3] provides evidence of the 
schematic phase diagram in Fig.QJa). 

Three types of Skyrmions can be continuously trans- 
formed into each other and the transition between them 
is identified as the second-order transition. The cylin- 
drical singularity of the component Vt^ and ^4. in the 
2D Skyrmion continuously deforms to sphere and ring- 
shaped ones and then the texture changes into the 3D 
Skyrmion. This is confirmed by the imaginary time evo- 
lution of the GP functional (JTJ) , which demonstrates that 
the 2D Skyrmion state becomes unstable in the vicinity 
of k z /k± ~ 1 toward the 3D Skyrmion. The 3D Skyrmion 
state also becomes unstable in the region of k z /k± > 1, 
where the singular line and sphere of described in 
Fig. EJJa). Hence, the 2D Skyrmion solution smoothly 
transforms to ID Skyrmion through the 3D Skyrmion 
as k z /k± varies from +0 to +00, and the 3D Skyrmion 
texture spontaneously appears at 

We emphasize that the smooth transition behavior be- 
tween Skyrmions ensures the stability region of the 3D 
Skyrmion against finite c\ ^ 0, which breaks the SU(2) 
symmetry in Eq. [TJ This is because the stability of ID 
and 2D Skyrmions in the limit k z — > and k± — > stays 
unchanged in Ci ^ 0, [3(| . In fact , within < Ci /cq < 1 , 
we numerically confirm that the winding number W3D in- 
creases as approaching k z /k±— >1 and the 3D Skyrmion 



^3D = 

lized. 



U^(r,2 K r)^ ] with (* f 
In the case of ci/c < 0, 



(0) ,t,(°) 



= (1,0) is stabi- 
the 3D Skyrmion with 



(#f *j UJ ) = (1, 1) is also stabilized. 

Skyrmions in other systems. — The stability of the 3D 
Skyrmion in other spinor systems is worth a mention in 
passing. For instance, we find that the 3D gauge field 
in Eq. (J3J) cannot stabilize the 3D Skyrmion in the polar 
phase of a hyperfinc spin F = 1 spinor BEC, which is the 
"knot" soliton. This is because the helical modulation of 
the polar OP is not degenerate for the 3D direction of 
the modulation vector In fact, the OP manifold 

in the polar phase reduces to an SO (2) symmetry in spin 
space. Hence, since the helical modulation vector along 
the direction of the £0(2) rotation axis cannot gain the 
single particle energy Eq, only the 2D or ID Skyrmion 
can be stable in this system. In the ferro mag netic phase 
of an F— 1 spinor BEC, the ID Skyrmion [20| and the 2D 
half-Skyrmion 37[ are proposed. However, the problem 
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on the stability of the 3D Skyrmion remains as nontrivial. 

Conclusions. — Here, we have demonstrated that the 
3D Skyrmion spontaneously appears as the ground state 
of two-component BECs coupled with a non-Abelian 
gauge field, which is a 3D analogue of the Rashba spin- 
orbit coupling. The appropriate gauge field and the 
spherical density distribution due to the trap potential 
are necessary to stabilize the Skyrmions. Upon squashing 
the 3D gauge field to ID or 2D shape, the 3D Skyrmion 
continuously undergoes a change to ID or 2D Skyrmion. 
We computed the ground state phase diagram with the 
winding numbers and helicity, which is covered by three 
types of Skyrmions. All of the emerged Skyrmions are 
understandable with the concept of the helical modula- 
tion of the order parameter. In addition to the SU(2) 
symmetric interaction, we also consider the stability of 
Skyrmions against an SU (2) nonsymmetric interaction. 
Then, we confirmed that the asymmetric term does not 
alter the phase diagram. However, the interaction in the 
presence of a synthetic gauge field depends on the detail 
of the scheme to imprint it. The complete phase dia- 
gram in the more realistic situation remains as a future 
problem. 
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Note added. — After we submitted this paper, we be- 
came aware of two papers [Z(], 41 1 which discuss the 
ground state properties under the same gauge field. 
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